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Three- wave interactions in plasmas are described, in the framework of kinetic theory, by the quadratic response 
tensor (QRT). The cold-plasma QRT is a common approximation for interactions between three fast waves. 
Here, the first-order thermal correction (FOTC) to the cold-plasma QRT is derived for interactions between 
three fast waves in a warm unmagnetized collisionless plasma, whose particles have an arbitrary isotropic 
distribution function. The FOTC to the cold-plasma QRT is shown to depend on the second moment of the 
distribution function, the phase speeds of the waves, and the interaction geometry. Previous calculations of 
the rate for second harmonic plasma emission (via Langmuir-wave coalescence) assume the cold-plasma QRT. 
The FOTC to the cold-plasma QRT is used here to calculate the FOTC to the second harmonic emission 
rate, and its importance is assessed in various physical situations. The FOTC significantly increases the rate 
when the ratio of the Langmuir phase speed to the electron thermal speed is less than about 3. 
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I. INTRODUCTION 

Plasma emission, which is the generation of radiation 
at multiples of the local electron plasma frequency f p , 
is observed in various locations and phenomena in our 
solar system; these include type Ilir— and type 
solar radio bursts in the solar corona and interplanetary 
medium, and terrestrial foreshock emission.^ Although 
several mechanisms have been proposed for these emis- 
sions, such as linear mode conversion^— and cyclotron 
maser emission^ they are generally attributed to three- 
wave interactions between Langmuir, transverse, and ion 
sound w&ves£12r— Three-wave interactions include the 
coalescence of two waves to give a product wave, and the 
decay of one wave into two product waves. These pro- 
cesses occur due to the nonlinear response of the plasma 
medium to the wave fields. In kinetic theory, the re- 
sponse of a plasma to an electromagnetic disturbance is 
described by linear and nonlinear response tensors.— ~— 
Assuming the plasma response to be weakly nonlinear, 
induced plasma properties, such as the induced charge 
and current densities, can be expanded in powers of the 
amplitude of the electromagnetic field: this is termed 
the "weak-turbulence expansion" 1 17 i 18 i 20 The quadratic 
response tensor (QRT), defined as the coefficient of the 
second-order term in this expansion, describes the re- 
sponse of the plasma to two fields beating simultaneously 
to produce a third wave; therefore, the QRT is the rele- 
vant response tensor for three-wave interactions. 

The general form of the QRT involves integrals over 
the velocity distribution function of the plasma parti- 
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cles. These integrals are difficult to evaluate exactly, and 
so the integrand is often approximated before performing 
the integrations, where the approximation that is made 
depends on the dispersion relations of the three wave 
modes involved. For interactions between three "fast" 
waves, that is, waves with a phase speed much greater 
than the thermal speed, the "cold-plasma approxima- 
tion" to the QRT is often made, in which thermal effects 
are neglected in the description of wave couplingj 17 i 18 
This approximation lends itself to a simple derivation 
of the response tensor, and the resulting expression can 
be readily used in practical applications. However, the 
cold-plasma approximation becomes less accurate as the 
phase speed of one or more of the interacting waves ap- 
proaches the thermal speed, and its range of validity is 
poorly defined. 

Conversely, the QRT has been evaluated exactly for a 
thermal plasma by Percival and Robinson—, in terms of 
generalized plasma dispersion functions. 2 — The resulting 
exact QRT accurately describes wave coupling in three- 
wave interactions between any wave modes. Although 
the use of the exact QRT is desirable, the expression 
is too cumbersome to apply analytically, and even its 
numerical evaluation presents difficulties as described in 
Sec. II. 

Due to the inaccuracy of the cold-plasma QRT at low 
phase speeds and the difficulty of applying the exact 
QRT, accurate approximations to the QRT are needed 
for a proper treatment of three-wave interactions. To 
this end, we derive here the first-order thermal correction 
(FOTC) to the cold-plasma QRT for a warm collisionless 
unmagnetized plasma, whose constituent particles have 
an arbitrary isotropic velocity distribution. The approx- 
imate response tensor derived in this paper, which is the 
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sum of the cold-plasma QRT and its FOTC, has the ad- 
vantages that it is more accurate than the cold-plasma 
QRT and more tractable than the exact QRT. This ap- 
proximation is valid for interactions between three fast 
waves, such as the Langmuir-wave coalescence process 
that generates second harmonic emission. We find that 
the FOTC depends on the second moment of the distri- 
bution function, the phase speeds of the waves, and the 
interaction geometry. 

First-order thermal corrections to the quadratic re- 
sponse have been derived before in the literature, but the 
cases treated are not suitable for modeling the Langmuir- 
wave coalescence process in space plasmas. For exam- 
ple, a FOTC to the cold-plasma longitudinal quadratic 
susceptibility was previously derived by Sitenko— ; how- 
ever, this quantity only describes interactions between 
three fast electrostatic waves and is thus inadequate 
for treating second harmonic plasma emission, in which 
electromagnetic transverse waves are produced. The 
FOTC to the cold-plasma QRT was derived in Ref. H 
for a Maxwellian velocity distribution of the plasma 
particles, but space plasmas are commonly observed to 
have power-law tails which must be modeled by a non- 
Maxwellian distribution, often the generalized Lorentzian 
(or "kappa") distribution! 24 i 25 The expression that we 
derive is equivalent to that in Ref. [23| when the velocity 
distribution is Maxwellian, but allows the treatment of 
three-wave processes in non-Maxwellian plasmas. 

The rate of a three-wave interaction is dependent on 
the strength of the coupling between the waves, which 
is described by the QRT. Until now, the rate of sec- 
ond harmonic plasma emission via Langmuir-wave coa- 
lescence has been calculated assuming the cold-plasma 
approximation.^— We use the FOTC to the cold- 
plasma QRT to calculate the FOTC to the rate of 
Langmuir-wave coalescence, and assess its contribution 
to the total interaction rate in various situations. There 
is a significant increase of the rate when the ratio of the 
Langmuir phase speed to the electron thermal speed is 
less than about 3. 

The paper is structured as follows. In Sec. II, the the- 
ory of response tensors and their derivation is described. 
In Sec. Ill, we derive and discuss the FOTC to the cold- 
plasma QRT. The FOTC to the rate of second harmonic 
emission is derived in Sec. IV, where its importance in 
various situations is also analyzed. 



II. THEORETICAL CONTEXT 

Nonlinear plasma response tensors are defined by ex- 
panding an induced plasma property, such as the in- 
duced current density, in powers of the amplitude of the 
Fourier transformed electromagnetic field; this is termed 
the weak-turbulence expansion. The QRT is the coef- 
ficient of the second-order term in this expansion. On 
choosing to describe the electromagnetic field by the vec- 
tor potential A in the temporal gauge, the induced cur- 



rent density is given in Fourier space by (e.g., Ref. [l8l ) 



(i) 



where 



4 n) (k) = I d\^a ijlh ... jn (k,k 1 ,...,k n ) 
x A h (ki)A h (k 2 )---A jn (k n ). 



(2) 



In Eqs. ([T]) and ([2]), k m collectively denotes uj m and k m 
for the roth wave, and is the nth-order convolution 

integral given by 



d 4 ki d 4 k 2 d 4 k n 



(2tt) 4 (2tt) 4 (2tt) 4 
x (2n) i 5 i (k -h- 



kn)i 



with 



and 



d k = efcjeTfc 



S 4 (k) = S(cu)S 3 (k). 



(3) 
(4) 
(5) 



The nonlinear response tensors are most commonly cal- 
culated via the Vlasov equation, which relates the dis- 
tribution function to the wave fields for a collisionless 
plasma. Solving the Vlasov equation by employing a 
weak-turbulence expansion of the distribution function 
and expressing the induced current as a moment of the 
distribution yields the general form of the QRT (e.g., 
Ref.[2l|), 



a ij i(k,k 1 ,k 2 ) 



2m 2 







fc l + (|fe| 2 - W 2 / C 2 )^ 



x [I(k)6jt + k 2 jJi(k,k 2 ) + kuJj(k,h) 
+ (fci • k 2 - LOiUj 2 /c 2 )Kji(k, ki,k 2 )] 

+ {i,k) O ft,**) + (*»*) O (l,k 2 )\. (6) 



where (i,k) ■<-> (J,ki) represents the additional terms 
generated from those written explicitly by interchang- 
ing k and k\ and the associated tensor indices, and the 
integrals I(k), Ji(kx,k%), and Kij(k,ki,k 2 ) are given by 



/(*)= / d 3 vf(v) 



1 



— k ■ v' 



Mh,k 2 ) = / d 6 vf{v) 



(wi - fci • v)(u)% - k 2 ■ v) 
K l3 (k,k u k 2 ) = Jd 3 vf(v) 



(lu — k ■ v)(lui — ki ■ v)(uj 2 — &2 • v) 



(7) 
, (8) 

(9) 
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It was noted in Ref. 

so that the partial derivative 
dKji/dki cannot be taken without ambiguity since one 
of the wave vectors in the integrand is no longer indepen- 
dent of the other two. The integral 



Liji — 



d 3 vf(v) 



ViVjVl 



(lo — k ■ v) 2 (uji — ki ■ v)(u 2 — k 2 ■ v) 

' (10) 



must therefore be evaluated directly. 

Thermal effects in the wave coupling are often ig- 
nored for interactions between three fast waves; neglect- 
ing these effects is known as the cold-plasma approxima- 
tion. The cold-plasma QRT may be calculated by substi- 
tuting f(v) = n5 3 (v) into Eqs. (|71)- (fT0)) then evaluating 
the integrals. The integrals J, K, and L vanish, and I(k) 
is replaced by n/uj. This leads to 



(cold) 

a) 



q 3 n 
2^? 



h5 3 i kij8. a k 2 i5j 



U) 2 



(11) 



The cold-plasma approximation is made in order to 
simplify the mathematical analysis, but the neglected 
thermal effects may become significant in some circum- 
stances. 

Percival has calculated the integrals in Eqs. (jTl)- (fl"U)) 
exactly for a plasma in which the particles have a 
Maxwellian velocity distribution. 2 - This analysis yields 
an expression for the QRT in terms of generalized plasma 
dispersion functions^ but calculating the interaction 
rate (as described below) is problematic because it in- 
volves integrals of the squared modulus of the QRT con- 
tracted with the relevant polarization tensors. Due to the 
large number of terms in the exact QRT and the possi- 
ble occurrence of catastrophic cancellation it is infeasible 
to use this response tensor directly in the calculation of 
rates. 

Once the response tensors are known, the emission 
and absorption of waves can be studied. These processes 
are often treated semiclassically j 16 ' 18 the waves are inter- 
preted as a collection of wave quanta with momentum hk 
and energy H\ujM(k)\. The occupation number Nm(k) is 
introduced, being defined as the number density of wave 
quanta within the elemental range d 3 k of k. This quan- 
tity is related to the energy density per unit volume of 
fe-space, Wm(&), by 



N M (k) = 



W M (k) 
haj M (k) ' 



(12) 



In a coalescence process the current j( 2 '(k) induced by 
the simultaneous response of the plasma to two wave 
fields Ap(fci) and Ag(k 2 ), given by Eq. ([2]), is the source 
of a third wave field A^,/ (&), where k = k\ + k 2 as im- 
plied by the delta function. Assuming the random phase 
approximation, the rate equation for the wave mode M 



dN M 
dt 



in the three-wave interaction -P(fci) + Q(k 2 ) —> M(k) is 
given by (e.g., Refs. and HI) 

f d 3 fcx f d 3 k 2 

x {N P (k 1 )N Q {k 2 ) - N M (k) 

x [N P (ki) + N Q (k 2 )]}- (13) 

Alternatively, the rate can be expressed in terms of Tm, 
the effective temperature for the wave mode M, as 

d 3 k 2 



dT M _ f d 3 k x 



rU M PQ(k, ki,k 2 ) 



(2tt) 3 J (2tt) 3 
x {Tp(k 1 )T Q (k 2 )-T M (k) 
x [T P (fc 1 )+r Q (fc 2 )]} 
x u3 M (k)/(fiujp(ki)u;Q(k 2 )), 



(14) 



where the effective temperature is related to the occupa- 
tion number by 

T M (k) = hw M (k)N M (k). (15) 

The equation for the interaction probability umpq is 
(e.g., Refs. [Hand [T 



UMPQ(k, h,k 2 ) 



Ah R M {k)R P [k 1 )R Q {k 2 ) 



e 3 \Ld M (k)ujp(k 1 )ujQ(k2)\ 



x \aMPQ(kM, kpi, kQ 2 )Y 



x (27r) 4 S i (k M - fcpi 



kQ 2 ), 



(16) 

where Rm is the ratio of electric to total energy in the 
wave mode M , and 

aMPQ(kM,kpi,k Q2 ) =aiji(k M ,kpi,kQ 2 ) 

x e* m (k)e P j(ki)eQi(k 2 ), 
(17) 

with eM(k) the polarization vector for the wave mode 
M. The quantity ku collectively denotes u)]\i(k) and fc, 
and similarly for kp\ and kQ 2 . In Eq. (|16l) . the delta 
function in wave vector is interpreted in the semiclassical 
description as expressing conservation of momentum, 



k = ki + k 2 , 



(18) 



and the delta function in frequency expresses conserva- 
tion of energy, 



u)M(k) = wp(fci) + uj Q (k 2 ), 
where the common factor of h is omitted. 



(19) 



III. FIRST-ORDER THERMAL CORRECTION TO THE 
COLD-PLASMA QUADRATIC RESPONSE TENSOR 

In this section, the general quadratic response ten- 
sor (QRT) in Eq. ^ is approximated by deriving the 
first-order thermal correction (FOTC) to the cold-plasma 
QRT given by Eq. {TTJ, We discuss the expression ob- 
tained for the FOTC, including its importance and va- 
lidity, for different interactions. 
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A. Derivation 



This gives 



First the resonant denominators in Eqs. (I7|- (|10p arc 
binomially expanded in powers of k m ■ v /u m (which may 
be expressed as v^/v^m where || is with respect to k m ) 
using the binomial expansion (1 — x)~ l = X)^Lo x " • ^ ne 
v$ m in the denominator of the expanded quantity en- 
sures that this expansion of the integrals is convergent, 
as shown below. We show that this expansion recov- 
ers the cold-plasma QRT plus additional terms of order 
(V/vc/,) 2n relative to the cold-plasma terms, where is 
the phase speed, V is the thermal speed of the particles, 
and n is an integer. The n — 1 terms are called the 
FOTC to the cold-plasma QRT. 



1(A) =- / d 3 vf(v) 



k ■ v ( k ■ v 



O 



k ■ v 



i~ J d 3 vf(v) + h j d 3 vf(v)v q 
^ d 3 vf(v)v r v s , 



Ji{k,ki) = / d a vf(v)vi 



k ■ v ki ■ v 



ui 



O 



k ■ v 



i— [d 3 vf(v)v i + ^-(^+ 1 ^ 
x / d 3 vf(v)viV s , 



(20) 



(21) 



(22) 



(23) 



Two additional assumptions are made to derive the 
FOTC: the first is that the distribution function is 
isotropic, and the second is that the thermal effects for 
each of the three wave fields can be treated equally, which 
will be discussed after the derivation. 



Firstly, we perform a binomial expansion of the reso- 
nant denominators of the integrals given by Eqs. (|7|)- (fT0l) . 



K ij (k,k 1 ,k 2 ) 



L ij i(k 1 ki,k 2 ) 



LOUJ\L02 



d 3 vf(v)viVj 



1 + 



' k ■ v 



- — / d 3 vf(v)viVj, 

UJU>lU) 2 J 



UJUJIUJ2 



d 3 vf(v)v l v j vi 



O 



k ■ v 



to 



UJIjJ\tjJl 



d 3 vf(v)viVjVi. 



(24) 
(25) 
(26) 



(27) 
(28) 



We let Vi — 5i q v q and Vi = Si r v r in the first and sec- 
ond integrals in Eq. (|2Uf respectively. Then substitut- 
ing the binomially approximated /, J, K, and L given 
by Eqs. (|2T))) - (j2"5)) respectively into the general QRT in 
Eq. ([6]) gives 



2m 2 



k l + {\k\ 2 ~Lo 2 /c 2 ) 



k r k s Sji k2jSir ( k s 

U> 3 UJUJ2 V oj 



dh 

W2 



1 d 6 vf{v) + H J — i H ^ AJv) 

\ oj z UJ0J2 (*JU>i ) 

(fel • k 2 - LUlL0 2 /c 2 ) 



k u 5 



it 1 k s k\ s . . 
— — h \B rs {v 



UJUJ1UJ2 



kiBfliv) + {\k\ 2 - lo 2 /c 2 ) C i:jl (v) + (*, k) o (j, h) + (i, k) o (I, k 2 ) 



(29) 
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Here wc have defined 



Ai(v) = / d 3 vf(v) 



(30) 
(31) 
(32) 



We now reexpress the integrals over velocity space in 
terms of moments of the distribution function f(v). The 
distribution function is defined such that the particle 
number density n is given by 



B ij( v ) = J d 3 vf(v)viVj, 

C ijl( v ) = / d 3 vf(v)ViVjVl. 



n = J d 3 vf(v). 



The moment of a quantity Q(v) is defined by 

(Q(t>)> = - [ d 3 vf(v)Q(v). 
n I 



(33) 



(34) 



We make the assumption that the distribution function 
is isotropic, hence f(v) = f(y) where v = (v 2 + v 2 + 
v 2 ) 1 ^ 2 and so / is an even function of v x , v y , and v z . In 
Cartesian coordinates, the tensor indices run over x, y, 
and z. The integrals Ai and Ciji then vanish because 
every choice of indices gives an integrand that is odd 
in one or all of the variables v x , v y , and v z , and the 
integration limits are symmetric about the origin. By 
symmetry in v x , V y , and v z all the diagonal components 
of Bij are equal, as are the off-diagonal components. The 
off-diagonal components of the Bij vanish due to oddness 
of the integrand, so B^ = BSij where B = B xx = B yy = 
B zz . To calculate B we note that B = (B xx + B yy + 
B zz )/3 and so 



Bij= 6 f J d 3 vf(v)(v 2 x +v 2 y +v 2 z ), 



n(v 2 )6i 



(35) 
(36) 



from Eq. ((Ml) . 

Substituting Eqs. ([33]) and ([36)1. and A q = Cyj = 0, 
into Eq. © and simplifying gives 



Otijl 










)- 




UJ 


M + 


hi 


w 2 / 





) 8ji_ kj (fci • k 2 - uj^/c 2 ) 
UILU1LU2 



+ (i, k) o (j, ki) + (i, k) O (I, k 2 ) 



(37) 



where the arguments of the response tensor have been 
omitted for brevity. On performing the derivative in 



Eq. ([3"7|. using d\k\ 2 /dki = 2fcj and dkj/dki = we 
have 



atiji 



q 3 n J kjSji 
2m 2 I uj 



\ k \ 2 8ji , k 2 j (h , k 2t 



+ 



kil ( kj k\j 



WWl \UJ OJi 



UJ° UJUJ 2 \U1 UJ 2 

„2' 



[ki ■ k 2 - CJ1W2/C ) hSj 



+ (\k\ 2 - UJ 2 /C 2 ) ( 2ki l jl J- k2]5a 4- 



+ (j, fc) o (j, fci) + (i, k) o (Z, fc 2 ) 
Rearrangement and factorization of Eq. 



W 2 W 2 W 2 Wl 



(38) 



g 3 n 



h S 



iOjl 



3 



Ifel 



2(|fc| 



yields 

w 2 /c 2 ) 



2m 2 ^ uj 

fel ■ fe 2 - WiW 2 /c 2N \ fci (5,7 



|fc| 2 - w 2 /c 2 /fc 1/( 5, 



+ 



k 2 j5ti 



UJ 



kukj 



k%jh 



U> 2 J UJ \ UJUJ\ CjJi UIUJ2 

+ (i, k) <-> (j, fci) + (i, fc) 4-> (Z, fc 2 ) 



k 2 jk 2 i 
oj 2 

(39) 



The first term inside the braces and its interchanges are 
identified as the cold-plasma QRT given by Eq. (fTTj) : the 
remaining terms are the FOTC to Eq. (fTTj) . denoted by 

Aa>iji. That is, 0^7 = a:^° ld ^ + Aa.yi. In the nonrel- 
ativistic (c — > 00) limit, the FOTC to the cold-plasma 
QRT, for an isotropic particle velocity distribution, is 



Aotijl = 



q 3 n (v 2 ) ( 


r/3|fe| 2 


2m 2 3 t 


A ^ 2 


/ kuSij 
x - 


k2jS t i 
\ — 


V w i 


UJ 2 


fcoi'fc/ 


k2jk 2 i\ 


LOUJ2 


<4 J 



k 1 -k 2 \ hS-j 



UJ\UJ2 J uj 

ki { k\ikj 

UJ \ U1UJ\ 



UJ A 

uj 2 



+ (i, fc) f+ (j, fci) + (t, fc) O (f, fc 2 ) L (40) 



B. Discussion 

The particle temperature T is related to the second 
moment of the distribution function by the definition 
kgT = m(i> 2 )/3, where ks is Boltzmann's constant. For 
a Maxwellian distribution 



/(«) 



( 27r )3/2y3 



-« 2 / 2 y 2 



(41) 



the thermal speed is given by V = ^kBT/m, whence 
(i> 2 )/3 is replaced by V 2 in Eq. (FTO!) ; this reproduces the 
expression for Ao^-; in Ref. |23| . 

The terms in the cold-plasma QRT and FOTC are of 
order l/v$ and V 2 /vi respectively, from Eqs. (fTTj) and 
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(l40|) . where v$ here represents the phase speeds of any 
of the three interacting waves. The ratio of these or- 
ders, V 2 /vi, therefore determines the significance of the 
FOTC: if it is much less than unity then the cold plasma 
approximation will be accurate, but as it approaches 
unity the cold plasma approximation begins to break 
down and the FOTC must be included for an accurate 
description of the plasma response. 

From Eqs. © and p0 | -p8 ]l . it follows that higher- 
order thermal corrections to the cold-plasma QRT involve 
integrals of the form 



(42) 



where i m are tensor indices running over x, y, and z. For 
odd n, each choice of indices will give an odd power in 
at least one of the variables v x , v y , or u z , hence in this 
case the integral will be zero. Therefore, the nth-order 
thermal correction will be O {(V/v^) 2 ™} relative to the 
cold-plasma terms. For the expansion of the general QRT 
by a binomial expansion of the resonant denominators to 
be convergent, one requires > V for each wave mode. 

Thermal effects from one or two of the participating 
waves may be more important than those from the other 
wave or waves. This is the situation for the Langmuir- 
wave coalescence process: although both Langmuir and 
transverse waves are fast, the thermal effects from the 
Langmuir waves will be more important since they have a 
significantly lower phase speed than the transverse wave, 
and so the terms involving the transverse wave phase 
speed may be neglected. Equation (14*01) may also be 
applied to the Raman scattering process L + T — s- T' 
which has been proposed for third and higher harmonic 
emission.— ~— However, the FOTC will not be as im- 
portant as in Langmuir-wave coalescence since Raman 
scattering involves two transverse waves and only one 
Langmuir wave. 



Langmuir waves L(k\) are assumed to be generated by 
an electron beam via a bump-on-tail instability, such that 
the phase speed of the Langmuir waves is approximately 
equal to the speed of the electron beam, i.e. v^i rj v^. 
Backscattered Langmuir waves L'(fc 2 ), with which the 
primary waves coalesce, are assumed to be the product 
of the decay process L —> L / + 5 , i 26 ' 28 ' 32 " — Since the mass 
of the ions is much greater than that of the electrons and 
the QRT has a m~ 2 dependence for each particle species, 
the ionic contribution to the QRT is neglected. 

We first outline the derivation of the cold-plasma in- 
teraction probability for Langmuir-wave coalescence (see 
e.g. Ref Il8l). Transverse waves have ex • kt = and 
hence e T ■ k t = 0. In the case of no spatial damping, kt 
is real and so e T -Kx = 0. On contracting the cold-plasma 
QRT in Eq. (TJT]) with e^ i (k T )e L j(k 1 )e L n(k2), i.e. with 
e^ ri K\jK2i, where k = fe/|fc|, we have 



(cold) / , , , s e n e ( G T ■ Ki e T ■ K2 i 

a TLL !(k T ,k 1 ,k 2 ) = 7^ ( — + — ); (43) 



2m 2 



V<j>2 



"01 



hence, using \x + y\ 2 = \x\ 2 + \y\ 2 + 2Re{x*y}, we find 



(cold) 
"-TLL' 



"4m4 
+ 2Re 



\e? T ■ k 2 \ 



u 4>2 u <t>i 
(e* T ■ Ki)*{e* T ■ k 2 ) 



(44) 



When the polarization of the transverse waves T is of no 
interest, an average over the two initial states of polar- 
ization and a sum over the two final states of polarization 
is performed. This leads to the replacement 



IV. FIRST-ORDER THERMAL CORRECTION TO THE 
SECOND HARMONIC EMISSION RATE 



\e T ■ ki\ 



1 I |2 
- \K T X Ki\ 



(45) 



In this section the FOTC to the cold-plasma QRT, 
given by Eq. (j40j, is applied to the rate calculation of 
second harmonic plasma emission via Langmuir-wave co- 
alescence. The ratio of the FOTC rate to the cold-plasma 
rate is derived, and its range of values is calculated for 
different physical situations to assess the importance of 
the FOTC. 



Derivation 



and similarly for 



that Rl 



\efp ■ K2\ 2 ■ Making the approximations 



and lot 



i 

2 ' 

as 



2ui p , and since Rt 



iu ill, ~ 5> L Pi T alLKA oiuuc -ii/r — 

we have the interaction probability for a cold plasma 



(cold) 
l TLL' 



4 £ 3 W 3 



(cold) 
TLL' 



(2 7 r) 4 ci 4 (A: T -fc 1 -fc 2 ). (46) 



The FOTC to quantities in this section will be de- 
noted by the prefix A, such that x rj x^ cold J + Ax for 
some quantity x. We define the electron thermal speed 
by v e = v/ (v 2 ) /3, where the angle brackets denote the 
moment of the electron distribution function. Primary 



To obtain AarLL', Eq. (|40|) is contracted with 
e Ti KijK2i- Grouping Aoltll' by order in v^t gives 



ActTLL' = Aa 



(0) 

TLL' 



Aa 



(i) 

TLL' 



Aa, 



(2) 

TLL' 



(47) 
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where 



(o) _' 

TLL' 



2ml 



3e T -K 2 2(e T • k 2 )(«i • « 2 ) 



(i) 

TLL' 



e 3 n e 



2(e T ■ Ki)(ki ■ n 2 ) 3e T ■ Ki 
2(k t ■ K 2 )(e T ■ Mi) 



(48) 



2m2 



(k t ■ K 2 )(e T ■ k 2 ) + (kt • ^i)(e^ • «i) 
2(k t • «i)(ej • K 2 )' 



A, 



* (2) 

TLL' 



e 3 n e 



2ml 



u <t>2 

(«i • K 2 )(e T 



«0i 

K 2 )(e T ■ k 2 )' 



From 



Aa 



(i) 

TLL' 



V<f, 2 

these equations, 
= O (v 2 1- 2 



(49) 



(50) 



Aa 



(o) 

TLL' 
(2) 



/v^v^t), and Aa ( TLL , 



= 0(v 2 /vl 
O (v 2 /v b v 2 T 



The electron beam speed is typically less than a few 
tenths of the speed of light and v^t > c, so to first order 

Aa T LL' ~ Aa^ LL ,. 

Including the first order thermal correction in the non- 
linear response implies cxtll' — &tll> + Ao:tll' ■ Hence, 
to the next order after the cold-plasma term, 



\OtTLL' 



a 



(cold) 
TLL' 



-2Re 



Aa. 



(0) 

TLL' 



(51) 



So, the second term on the right hand side of Eq. (ISTj) is 
the first order correction to \ctTLL'\ 2 , which is then 



A(\a TLL ,\ 2 ) =- 



2ml 



Re 



3\e T ■ k 2 \ 2 [2{e* T ■ k 2 )(ki ■ k 2 ) + e* T ■ k x ] {e* T ■ k 2 )* + 3(e T ■ k 2 ) (e* T ■ KiY 



[2(e T ■ k 2 )(ki ■ k 2 ) + e* T ■ Kjj (e T ■ «i) 



+ («i,W0i) (n 2 ,v,p 2 ) 



(52) 



From Eqs. <JT6j), ((SU, and ((52]), the FOTC to the inter- 
action probability is 

Autll' « -^A(\a T LL'\ 2 ) (2Tr) 4 S 4 (k T -k 1 -k 2 ). 

(53) 

The FOTC to the interaction rate is next calculated in 
terms of the effective temperature using a modified sys- 
tem of spherical coordinates (as in Ref. |28| ). where the 
Langmuir wave numbers k\ t2 take on positive and nega- 
tive values, and the polar angle 9 ranges from to ir/2. 
The angle between the primary and backscattered Lang- 
muir waves is assumed to be greater than 7r/2, hence 



sgn[k 2 ] = — sgn[fci]. If the back-reaction T — > L + L' is 
ignored, the FOTC to the rate of Langmuir-wave coales- 
cence is given approximately by 

A <gr f . )= _2_ f d 3 k! f d 3 k 2 
dt [ T) frwpj {2ir) 3 J (2tt) 3 
X AuTLL'{k T , fcl, k 2 ) 

x T L {k 1 )T L ,{k 2 ). (54) 

We assume that u>i 2 w oj p , whence v^i 2 = u> p /ki :2 , to 
simplify the integrand. The delta function S 3 (kx — k% — 
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k 2 ) is used to integrate over d 3 k 2 . Thus on integration, 



k-2 



k T - fei, 



e T ■ n 2 



e T ■ Ki, 



k 2 (ki,k T 
kr cos ip — k\ 
k 2 (k 1 ,k T ) 



(55) 
(56) 

(57) 



where ip is the angle between the T and L wave vectors, 
and 

k 2 {k 1 ,k T ) = -sgn[fei] (k\ + k% - 2k x k T cos ip) 1/2 . (58) 
The rate then becomes 



with 

g{ki,k T ) 



1 V)z 2 t ,m 



2 L, 



d 3 k 1 g(k 1 ,k T ) 



x S [io T {k T ) - Wi(fei) - wl(A; t - fei) 
xTLik^TvikT-k^le^-ml 2 , 



(59) 



3fc° 



[fe(fcl,feT)] 2 
+ fc2[fc 2 ( fclj fc T )]2 

fcf 



2fc? 



- 2 



[fe(fci,fe T )] 2 
(fcr cos ip — ki) — 2ki 



k\ (kx cos ip — k\) 
[fe(fci,fc T )] 2 
4fci (fey cos ^ — fci) 

2fc 1 [fc 2 (fc 1 ,fc T )] 3 



k 2 (k u k T ) 



3[fc 2 (fci,fe T )] 4 . 

(60) 



Expanding Eq. (|60l) , simplifying, and factorizing leads to 



g(ki,k T ) = 



|fc 2 -[fc 2 (fci,fcr)] 2 } 2 

[fc 2 (fci,feT)] 2 
x {3fc 2 + 3[fc 2 (fc 1 ,fe T )] 2 +2fc 1 

x (k T cos V' — fci)} (61) 
k\{2k\ costp — &t) 2 

[fc 2 (fel,fcT)] 2 

x (4fc 2 + 3fc£ - \k x k T cos-0). (62) 



On summing the states of polarization of the T waves 
in the final state and averaging over the initial states of 
polarization, • «i| 2 is replaced by \kt x k 1 | 2 /2 = 
(sin 2 ip)/2 in Eq. (J5T 



To simplify the delta function, the assumption 



Wi(fc) « w p + 3k 2 vl/2uj p 



(63) 



is made for both Langmuir waves, which is valid for 
k -C A^ 1 = (ve/uip)^ 1 . The dispersion relation for Lang- 
muir waves in a generalized-Lorentzian plasma (i.e., one 
in which the electrons have a kappa distribution) is also 
given by Eq. (|63|) in the limit k -C A^ 1 . 35 So, on sub- 
stituting Eq. (|6"3")l into the delta function in Eq. ([5U]) . we 
find 



£>[w T (fcT) - Wi(fel) - Wl(&T - fel)] 

f 3w 2 
= 5 < [ui T (k T ) - 2uj p ] - (2k\ + k% - 2k x k T cos% 

(64) 

Using 5(ax) = <5(x)/|a| gives 

S[u> T (k T ) - wi(fci) - uj L (k T - fei)] 



3w 2 |fci|/c2 



8 1 cos tp 



1 



2fcifc 2 



|^f[w T (fc T ) - 2w p ] 



(65) 



Evaluating the integrals over cos# and <p leads to the 
rate equation 



dk 1 g{ki,k T ) 
cos V>(fci, kx) ^1 + 

1 /? 

/3 sin x [l - cos 2 ?/>(fci , k T )] 



487re '7i 2 w 3 



x exp < p3 cos x 



x I n 



\k 2 (ki,k T )\J k 2 {kx,kT) 



\k 2 {k u k T )\ 



TLik^TMkukr)] 



(66) 



I 

where Iq is a modified Bessel function, x is the angle between fcx and the fey axis, and cos?/' satisfies 

COSf/^fcl, fc T ) 



_L_^ ; + 4_|| MM _ 2 .j) 
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(67) 



B. Discussion 



Substituting Eq. ([57]) into Eq. ||5JJ| yields 
|fci|fc T [2fci coaip(ki,kr) — k T ] 2 



g{ki,k T ) 



Mk^kr)} 2 
x [1 — cos 2 ip(k\, fcr)]- 



2w„ 



The cold-plasma interaction rate derived by Li ef 
is given by Eq. (|66[) on replacing g{k\, for) with 



(68) 

„7 28 



h(k\, kj 



uj 2 \ki\k T [2kiCOSip(ki,k T ) 
: 2^| Mkukr)] 2 
x [1 — cos 2 ip{ki, kr)]- 



The ratio R = g(k\, kr)/h(ki, kr) is then given by 



R(hr) 



2k T v 2 8 [u T (kr) - 2cj p 



3w c 



(69) 



(70) 



Since R is not a function of ki , Eq. ([TO]) gives the ratio of 
the FOTC to the cold-plasma interaction rate, indepen- 
dent of the integral over ki . On subsituting the minimum 



transverse wave number kro 
ratio can be expressed as 



R(k 7 



>V3/c into Eq. ((TTJ), the 



(—} ( — 

\ c J \k T o 



'TO 



(71) 



To obtain the emission rate, a particular Langmuir 
spectrum Iz,(fc) must first be assumed, after which the 
integral in Eq. (|66|) can be evaluated. For a Gaussian 
Langmuir spectrum, 



= cxp 



K 2 



-exp 



(fc - fc fc ) ; 



(72) 



the emission rate peaks at a transverse wave number 

&Tmax = &T0(1 + e) (73) 

where e = (A: 2 +k 2 )jk 2 D with fco = A^ 1 . Since the L' 
waves are produced by the electrostatic decay process 
L—>L' + S,kb~kf — fco ; 26 ' 27 where fc/ = uj p /vb and 
fco = 2wpD. 5 /3'u 2 with u s the ion acoustic speed. This 
leads to 



2^ s ^b 



2^| 
9v 4 



(74) 



in Eq. ([721). 



The FOTC to the cold-plasma rate of second harmonic 
emission has been derived by applying the FOTC to 
the cold-plasma QRT. This derivation is valid for both 
Maxwellian and generalized Lorentzian distributions of 
electrons; this is due to the FOTC to the cold-plasma 
QRT being valid for arbitrary isotropic velocity distri- 
butions, and to Langmuir waves having the same disper- 
sion relation for both distributions in the long wavelength 
{kX D < 1) limit. The resulting ratio R{k T ) of the FOTC 
to the cold-plasma emission rate, given by Eq. (|7Tj) . does 
not depend on the integration over k\ , so it is the same 
for all Langmuir wave spectra. 

We define the quantity R max to be the ratio R(kx) 
evaluated at kr = fcTmax in Eq. ([7T1) . We choose i? max 
to quantify the importance of the FOTC to the second 
harmonic emission rate. For most applications v e /c <C 1, 
and so the first term in Eq. (fTTj) is small. Importantly, 
Eqs. ([73"]) and ([74]) then imply that i? max depends mainly 
on the ratios Vb/v e and v s /v e , not on the individual 
speeds. As Vb/v e decreases, fcTma X increases, and hence 
-Rmax increases. Since v s /v e <C 1 unless T; ^ T e , the 
final term on the right hand side of Eq. ([74]) can be ne- 
glected, so R m ax decreases with increasing v s /v e . The 
dependence of i? max on Vb/v e is stronger than on v s /v e , 
which can be seen in Fig. [1] 

The significance of the FOTC to the cold-plasma rate 
of second harmonic emission is now assessed for differ- 
ent applications. In coronal type III solar radio bursts, 
typical parameters are Vb ss (0.2 — 0.5)c, v e ss 0.02c and 
v s = 1.5 x 10" 4 c. This gives k Tma Jk T0 = 1-02 - 1.003 
from Eqs. ([73]) and (71]), and hence i? max = 0.08 - 0.01 
from Eq. ([7T]) . However, Dulk et al— determined much 
lower electron beam speeds from their observations, rang- 
ing from Vb — (0.07 — 0.25)c, with an average of 0.14c. 
These values lead to the range fcTmax Ato = 1-16 — 1.04 
and i? max = 0.66 — 0.05, with an average of fcTmax/fcro 
1 .04 and i? m ax = 0.16. Thus, the second harmonic emis- 
sion rate may be well approximated by assuming a cold 
plasma for faster beams, but the FOTC becomes impor- 
tant for the slower electron beams measured by Dulk et 
al 

The electron beams responsible for significant radio 
emission are typically much slower in the "foreshock" 
regions upstream of shocks. Examples are Earth's fore- 
shock radio emissions, produced upstream of Earth's bow 
shock, and type II solar radio bursts associated with 
traveling shocks. At Earth's foreshock, n e «3x 10 _3 c, 
v s ^3x 10- 4 c, and v b /v e « 2- 10 (Ref.Hl), which gives 
fcTmax/fcro = 1.44-1.01 and i? max = 1.82-0.04. Knock 
et al— calculated, for interplanetary type II bursts, a 
maximum in the emissivity of second harmonic radiation 
where Vb/v e ss 3.5 for a thermal speed v e — 0.005c. In 
this case, taking « s = 1.3 X 10 _4 c leads to fcrmax/^ro = 



1.15 and R n 



0.63. These values of i? max indicate that 



the FOTC may be a significant contribution to the total 
rate in foreshock emissions, and can even exceed the cold- 
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0.1 0.2 0.3 

v /v 

s e 



FIG. 1. i? m ax versus (a) Vb/v e , where v a /v e = 0.1; and (b) 
v s /v e , where Vb/v e = 4. For both (a) and (b), v e = 0.003c. 

plasma contribution. Figure [5] shows the emission rate 
versus kr for typical coronal type III burst and Earth's 
foreshock parameters. Notably, the peak wavenumber 
^Tmax stays almost constant when the FOTC is added to 
the emission rate. 

Assuming the first term in R(kx) in Eq. (|7Tj) to be 
negligible, and also that v s /v e < 1, we obtain R > 1 for 
Vb/v e < 2.9. Thus, for sufficiently slow electron beams, 
the contribution from the FOTC exceeds the cold-plasma 
contribution to the emission rate. However, the assump- 
tion made in Eq. (|63|) that k^ <C A^ 1 , which corresponds 
approximately to Vb /v e S> 1 , is not satisfied very well for 
these slow foreshock parameters. Thus for small Vb/v e 
the expression for cos V^&i, &t) given in Eq. (|67|) . and 
hence the rate in Eq. ([66]) . will be less accurate. 

V. SUMMARY AND CONCLUSION 

Both the cold-plasma and the exact quadratic response 
tensor (QRT) describe three-wave interactions in which 
each wave is fast (that is, its phase speed is greater than 
the thermal speed). However, neither is ideal for the 
calculation of interaction rates: the cold-plasma QRT is 
readily calculable, but is only accurate where all phase 
speeds are much greater than the thermal speed; con- 
versely, the exact QRT provides an accurate description 




1.005 1.01 1.015 1.02 1.025 1.03 

Wo 




1.05 1.1 1.15 1.2 1.25 1.3 1.35 



V k T0 



FIG. 2. Second harmonic emission rate for (a) coronal type III 
burst parameters: v e — 0.02c, Vb/v e = 10, and v s /v e = 0.03; 
and (b) Earth's foreshock parameters: v e = 0.003c, Vb/v e — 3, 
and v„/v e = 0.1. Dashed lines are for a cold plasma while solid 
lines include the FOTC. 



of three- wave interactions between any waves, but its di- 
rect application to the calculation of rates is infeasible. 
The approximate QRT that we have derived here, which 
is the sum of the cold-plasma QRT and its first-order 
thermal correction (FOTC), overcomes these disadvan- 
tages since it is more accurate than the cold-plasma QRT 
alone but still permits a calculation of the interaction 
rate. It is also valid for arbitrary isotropic velocity dis- 
tributions. This approximate QRT is therefore suitable 
for modeling three-wave interactions in space plasmas, 
in which thermal effects are important for the interact- 
ing waves, and the velocity distributions are commonly 
non-Maxwellian. 

The rate of second harmonic plasma emission via 
Langmuir-wave coalescence has previously been treated 
with the cold-plasma QRT. Therefore, the resulting ex- 
pression is inaccurate where the phase speed of one or 
more of the waves is similar to the thermal speed. Using 
our result for the approximate QRT, we have derived the 
FOTC to the rate of second harmonic plasma emission. 
The ratio of the FOTC rate to the cold-plasma rate is 
easily calculated using Eq. (|7T1) ; it is only a function of 
the transverse wave number kx, and does not require an 
integral over the Langmuir wave number. The impor- 
tance of the FOTC to the emission rate is determined by 
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the ratios Vb/v e and v s /v e : the FOTC to the emission 
rate becomes larger compared to the cold-plasma emis- 
sion rate as both Vb/v e and v s /v e decrease. The FOTC 
to the cold-plasma emission rate is therefore important 
in foreshock emission, where the electron beam speed is 
not much larger than the electron thermal speed (within 
a factor of <~ 2 — 10). In the case where v s /v e <C 1, the 
FOTC to the cold-plasma emission rate is greater than 
the cold-plasma emission rate for Vb/v e < 3. 

Future work will involve deriving more accurate ex- 
pressions for the rates of the processes L — > V + S and 
L + S — > T involved in plasma emission. 
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